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ABSTRACT 

Two simple microwave radar techniques have Iven advanced that are potentially capable of pro- 
viding routine satellite measurements of the directional spectrum of ocean waves. One technique, the 
short-pulse" technique, makes use of very short pulses to resolve ocean surlace wave contrast fea- 
tures in the range direction; the other technique, the “two- frequency correlation" technique makes 
use of coherency in the transmitted waveform to detect the large ocean wavi contrast modulation as 
a beat 01 mixing frequency in the power backscattered at two closely separated microwave frequen- 
cies. A frequency-domain analysis of the short-pulse and two-frequency systems shows that the two 
measurement systems are essentially "duals", they eaeli operate on the generalized (three-frequency) 
fourth-order statistical moment of the surlace transfer (unction in different, but symmetrical ways, 
and they both measure the same directional contrast modulation spectrum. A three-dimensional 
physical optics solution tor the fourth-order moment is obtained for backscatter in the near vertical, 
specular regime, assuming Gaussian surface statistics. I he modulation spectrum is found to be given 
by the two-dimensional Fourier transform of the product of the joint surface height characteristic 
function and joint specular surface slope probability density function. Linearization in terms of the 
surface height and slope covariance functions yields a modulation spectrum that is directly propor- 
tional to the large wave directional slope spectrum evaluated in the direction of radar azimuth. A 
sample calculation with a model one-dimensional surface spectrum is carried out to indicate the ex- 
tent ol harmonic distortion ol the slope spectrum due to second-order nonlinear terms. For incidence 
angles 0 = 10-1 5°. the distortion is less than U) percent over the range of roughness conditions of 
practical interest. In principle, this distortion can be removed by iterative deconvolution. 
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DIRECTION Al SPECTRA Ol OCEAN WAVES FROM MICROWAVE 
BACKSCATTFR: A PHYSICAL OPTICS SOLUTION WITH 
APPLIC A 1 ION TO THE SHORT-PULSE AND TWO-FREQUENCY 
MEASUREMENT TECHNIQUES 

INTRODUCTION 

Two simple microwave radar techniques have been advanced that are potentially capable of 
providing satellite measurements of the directional spectrum of ocean wind waves on a routine, 
global basis. These two techniques, the “short-pulse” technique and the “two-frequency correla- 
tion" technique arc very similar in nature. Both approaches involve the detection of the modu- 
lation of ba k scattered power at oblique incidence angles caused by the modulation of "small 
scale scatterers” by the large underlying gravity waves. 

The small-scale scatterers-depending on the radar beam incidence angle— may be specularly 
reflecting wave facets, or Bragg-resonant diffraction elements; the modulation mechanism may he 
purely geometrical, js in the tilting of the small-scale scattered power pattern by the slopes of the 
large waves, or the modulation mechanism may be hydrodynamic in nature, involving the modi- 
fication of the small-scale waves by the atmospheric and large ocean wave flow fields (Alpers and 
llasseltnann, 1978). In the short-pulse technique, very short, wide-hand pulses are used to resolve 
the contrast modulation in range; the modulation spectrum is obtained simply by analog spectrum 
analysis of the envelope-detected backscatteicd signal. The two-frequency technique makes use 
ol coherency in the transmitted waveform; the large-wave contrast modulation is detected in the 
beat, or mixing frequency, of backscattered power at two closely separated microwave frequen- 
cies. flic short-pulse approach to a possible satellite ocean wave sensor was first suggested in a 
note by I'omiyasu (1971 ). The two-frequency approach, first proposed by Ruck, et al., (1972). 
has been extensively analysed by Alpers and Hasselmann (1978) who found it to be a feasible 
technique for aircraft and satellite implementation. Unfortunately no such similar analysis of the 
short-pulse technique exists in the literature, and as yet no “trade-off analysis of the two tech- 
niques has been made, although a beginning was made by Jackson (1974b). 
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As shown by Jackson (1 474b). ami as will be shown in the present work, the two technk|ues 
measure the same tiling namely, the directional contrast modulation spectrum of the large-scale 
ocean waves hi the direction of radar azimuth. High directional resolution is achieved through the 
lateral averaging of wave contrast features across the antenna beam spot on the surface, a distance 
which might be on the order of I Okm for a satellite radar system. The effect of this lateral aver- 
aging. when combined with spectrum analysis of the returned signal, is easily seen to In* that of a 
two-dimensional spectral analysis for contrast waves propagating in (or contrary to) the diicction 
ol radar look. In a nutshell, only those plane Fourier contrast waves that are travelling m the 
direction of radar look can survive the lateral averaging (Fig. lb). 

Two important developments have enhanced the potential for space application ot the two 
“directional wave spectrometers.” The first development (Jackson, I ‘>74) was the realization that 
the modulation signal spectrum to the residual background, random clutter* noise ratio (SNK. tor 
short) could be greatly improved in the two-frequency system by using wide-hand coherent wave- 
forms as opposed to the quasi- monochromatic waveforms originally figured m this technique. The 
second important development (Alpers and Hassclmann. I‘>78) was to show how Doppler filtering 
could be used on fast-moving platforms, also, to substantially increase the measurement SNR. 

In the present work, we will touch only briefly on general system considerations (e.g.. a rough 
analysis of the measurement SNR for the two techniques is carried out in the last section of the 
paper). I'lte primary aim of this paper is to present a solution for the directional modulation spec- 
trum m the restricted case of near-vertical incidence (t) < 10-20°). i.e.. in the case* of quasi- 
spcculai backscatter. Apart from the important limitation inherent in an assumption of Gaussian 
surface statistics, the physical optics solution presented is felt to be a rather realistic solution to 
the problem. The reason for putting faith in the solution is that in the specular sea backscatter 


•The expression "Rayleigh clutter" is used for the random, fluctuating pari of the backscattercd signal that arises 
solely as a consequence of waveform coherency. The superposition of monochromatic waves scattered from a 
larger number of independent scattering elements landomly distributed in range results in Gaussian statistics for 
the field (Central I unit Theorem) and classical Rayleigh statistics foi the envelope. In analogy with the conven- 
tional use of the word “clutter" for sea-echo interference with hard target detection, we use the word for the 
random fluctuations (hat interfere with the measurement of the modulation on the sea echo signal itself. 
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Finnic la. Overall Geometry and Coordinate Definition 



Figure lb. Illustrating Directional Resolving Power (In the 
Case Where Wave-front Curvature Is Not the 1 uniting Factor) 



regime, the sensitivity of the return to very small stale waves (e.g., capillary waves lor 2-3 cm 
radars) is small compared to the sensitivity to these waves in the L*ragg resonant diffraction regime 
at larger incidence angles <0> 30°). where the sensitivity is essentially 100'? (Wright. 1968). In 
the quasi-specula r backscattcr regime (Barrick, 1974). the entire sea spectrum up to the small 
diffraction-limited wavelets is contributing to the population of specularly reflecting facets, and 
there is, therefore, some reason to suppose that the small scale wavea-whose statistics via-a-vis the 
large wave profile are presently virtually unknown— will not enter into the problem in a too-critical 
way. The successful prediction of the average radar impulse response of the sea at vertical in- 
cidence by Jackson (1979) supports this supposition: in that work, assuming a geometrical optics 
model of scattering it was found that the (3cm) data could be explained by the non-Gaussian 
stati 'tics of the entire free gravity-wave ensemble without special consideration of the small waves. 
Thus, in the solutior we present here (a physical optics solution in the high-frequency limit), it is 
probably the Gaussian assumption for the large wave field that is the weak point. 

The above is not to say that larger incidence angles are to be avoided, the point is that theory 
for the hydrodynamic interactions that are the primary agents of the modulation at large angles is 
on very weak foundations (Wright, et al., 1978). A reliable model for large-angle backscattcr 
would have to be constructed on the basis of empirical modulation functions: presently there are 
insufficient data to define these functions (or functionals) over a range of sea states. 

In the following section we present a straightforward frequency domain analysis of the short- 
P’.ilse and two-frequency systems that allows lor arbitrary transmitted waveforms. The two sys- 
tems can be viewed as alternative systems lor detecting a modulation on noisy backscattcrcd sig- 
nals, a viewpoint inspired by the work of Par/en and Shircn (1956). It is shown that each system 
operates on the generalized three -frequency fourth-order statistical moment of the surface scatter- 
ing transfer function, M(k, k. Ak), in different but symmertical ways. M(k, k, Ak ) is ealeui; ted 
according to physical optics in three-dimensions, in the plane wave approximation, assumi ig a 
Gaussian sea. It is shown that M(k. k. Ak) - M , ( k . k) + M ( (k. Ak I and that the short-pulse* 
system modulation spectrum is the non-DC part of M , < k. *c>. and similarly for the two-frequency 
system with K-*Ak. In the limit k-*°°. M(k) has the form ot the two-dimensional I nurier 
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transform of the product of the joint surface height characteristic function and the joint specular 
surface slope probability density function. Linearization of the integral M, in terms of the height 
and slope covariance functions results in direct proportionality of to the directional large wave 
slope spectrum (a result anticipated). Harmonic distortion or smearing of the slope spectrum due 
to nonlinearity in can be accounted for accurately, and in principle rectified by an iterative 
deconvolution procedure, using the second-order expansion of M,. A sample (two-dimensional) 
calculation with a model wave spectrum is given to indicate the possible extent of harmonic dis- 
tortion. For angles 10° <0 < 15°, the distortion is < 30% provided a reasonably sufficient amount 
of roughness (i.e., total radar-effective rms surface slope). A brief analysis of the measurement 
signal-to-noisc is carried out in the last section of the paper, an analysis which includes the 
Doppler filtering concept of Alpers and Hasselmann (1978). The analysis indicates the feasibility 
of satellite measurements of directional wave spectra by either short-pu’se or two-frequency 
techniques. 


SHORT-PULSE AND TWO-FREQUENCY TECHNIQUES: THE GENERALIZED FOURTH 
ORDER MOMENT M(k, k, Ak) 


We consider the backseat ter of ;>n arbitrary (finite-energy) incident electric field waveform 
E 0 (t) whose Fourier transform is E 0 (k) where k = n/c is the propagation constant, v and c being 
the radian frequency and speed of light respectively. E„(k) is understood to be evaluated at the 
beam spot center x = (x, y) = (0. 0) on the mean surface z = 0. The surface transfer function for 
backscatter, Stk), is the ratio of the backscattered field harmonic component E s (k) to the incident 
field harmonic component, viz. 


S(k) 


E s (k) 

F. 0 (k) 


( 1 ) 


where the same holds in terms of frequency, v = kc. In principle, S(k) is obtained as the solution 
to the complex boundary value problem lor a unit amplitude incident monochromatic wave. If 
the scattering surface /. = f(x) is assumed to be homogeneous random process in the horizontal 
coordinate x, then S(k) is a complex random variable (approximately Gaussian, if the illuminated 
area is large compared to the scale of the spatial surface correlation, namely, the dominant ocean 
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wavelength) that obtain* lor an> rcali/ation of ftx). (In this work, we will not model the “slow- 
time.” Doppler evolution of S<kl; hut we will consider this aspect of the problem when we discuss 
the measurement SNR problem.) 


The short pulse detection system, diagrammed in figure 2a. consists of a square-law envelope 
detector and a spectrum analyzer. Hie ensemble average output of the spectrum analyzer can 'v 
written as: 


;j)> = (2n)* J I II, 


<0<uij» - (2ffi 3 / I ll,(u>)l* <1 l’(u>)l *> ilu; 
where II, is the Hh bandpass filter function and where 
<IP(cj)l 2 


l : > = J(S(v)S*(v - u>)S*(»'')S(»'' - u>)> • I „(f)l- 0 *(f - u>)l o*(*')l 0 ( 1 ’’ " wldrdr*'. 


The two-frequency detection system, diagrammed in Figure 2b. is a wide-band version of 
the original narrow-band (monochromatic) system first proposed by Ruck et al.. ( l*>72) and con- 
sidered by Alpers and llasselmann (l l, 7K). It consists of a two channel receiver and a cross- 
correlator. II, and II, are bandpass filters with center frequencies r, and n, separated by the 
modulation frequency An,, = i' ( - r,. I he K— filters are high-pass, l)( -blocking filters. The en- 
semble average output of the correlator can be written as 


<Q 12 > = (2 n)*J K, K 2 *<P, P 2 *>dw 

<P| Pj *> = ( 2n)*J" - u>)S {v )S(r»* - w)> • II, (i')ll, *(«' - w)H, *(n')ll 2 (r'' - w) (3) 


• I „(»')! „ *(*’- cell 0 *ti'')l 0 (n' - uMdndn'. 


The generalized fourth-order moment ot the transfer function tliat both systems operate 
upon will lie denoted by 

Mtk k. Ak) <Stk)S*(k - K)S*(k’)S(k' -k)> (4) 

where again we are using wavenumber k and (requeue) r interchangeably. I lie difference or 
mixing wave number is defined by 


n 
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Figure 2a. Short-Pulse System Diagram 



Figuie 2b. Two-Frequency System Diagram 
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Ak • k - k' 


It the scattering process d.e.. the surface impulse response) is described by a simple, phys- 
ically plausible mathematical model, to wit. a weakly modulated white (iaussian noise process, the 
functioning of the two systems for different transmitted waveform characteristics in detecting the 
modulation spectrum is easily analyzed (Par/en and Shiren, l l >5(>). We will hold off on this analy- 
sis until we have obtained a more meaningful physical scattering solution for Mtk. k , Ak). 

A THR1 l -DIMI NSIONAl PHYSICAL OPTICS SOLUTION I OR Mtk. k. Ak) 

A good solution* tor S(k) in the near-vertical, specular backscatter regime can be had with 
scalar physical optics using the so-called kirchhoff or tangent plane boundary values of the field 
in the physical optics integral The approach that we take in computing the fourth-order moment 
parallels the approach of Beckmann and Spi/zichino (|4(>3) to computing the second-order mo- 
ment <lS(k)l J >* That is. we will go to the lug*, frequency limit k-*°° after interchanging expecta- 
tion and integral operations. The alternative approach is to make the stationary phase* approxi- 
mation to the physical optics integral bet ore ensemble averaging (see. e.g.. Weissntan. 1474). 

We work under the simplifying assumptions ol (a) l ar /one, plane wave scattering. This 
assumption is good if we stay away from vertical incidence and keep to high altitudes: The phase 
front should be tlat over the height extent of the surface and over the azimuthal large-wave cor- 
relation distance, (b) Perfect conductivity. This is .1 tiiviul simplifying assumption (a Fresnel re- 
flectivity at normal incidence can be applied at the end), (c) Deep phase modulation, ko»l 
where o is the rms surface height, fins is entirely valid for microwaves ol 2-3 cm wavelength in 
any sea state, (d) Moderate bandwidth*. A high bandwidth in the short pulse system would be 
1% (Le., 2cm/2m). (e) Ciaussi in surface statistics. This is probably a poor assumption. The work 
of Jackson (147K) demonstrates the importance of the non-linear wave statistics in the average 
impulse response in the near-vertical. 


*1 xactly how good is this solution? A concrete, quantitative answer diould he possible by working on the wor k of 
Axline and lung ( 1 478) (numerical. Monic-l'arlo solutions ol integral equations) Brown ( l‘> 7 X) (composite, 
large and small scale surface method): and Jackson ( |474a) (high-frequency . parabolic correction to the 
kirchhoff tangent plane boundary values in the Stratton-Otu integral). 
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Let z a f(\ > describe the landoni surface. Following Wciwnan (1974), we write the physical 


optics integral, as 


Stk) a 


C k /cc i i e ukco ^f^ ) e -i2k«n9x d a 


<M 


where x. z form the plane of incidence 
0 is the angle of incidence 
S ■ (x, y) and d\ * dxdy 

* ±ik(2>r Rcosflr'exptikK) where R is the range to the beam spot center, \ * Q 
G( x ) = antenna power pattern, one way, G(Q) * I. 


Assuming interchangeability ol expectation and integration operations, the fourth order mo- 
ment (4) can tv written as the four-fold integral. 



G 4 ( e* *’ co *^* h e ‘ 2 9 k ’ * d x , 


• • d&a 


(b) 


In this, we use the short hand notation 


X 2 (Xj) * (x, . X,, Xj, x 4 ) 

k s (kj) « (k. - k ♦ k. - k / , k' - k) 

T -s lfi> * f<Xj). £<&,). s(x 4 ) 

The approximation to the product of the C^’s, folh>ws from the moderate bandwidth ay 

sumption, (d). The dot notation is the conventional vector inner product. 


The expectation 

0<2cosdk) * < e i2cos0k-f) 

is the characteristic function of the random loui-vector of surface heights f. If f normally div 
tributed with zero mean, then (I'isz, |9t>3)0 has the form 

0(2cosdk) * cxp(-2cos 2 0 ^ Rjjkjkj) (7) 


o 




when* the covariance matrix 


K ( j - a Ktjj-ij) (8) 

is a symmetric matrix and a symmetric function in each of the lag vectors At zero lag 

Rfj(Q) * (f 2 > * oJ <d) 

where u is the rms height of the surface. Let 

U = (u x . u. i » x ^ - x | 

v 5 <V X . V y ) = ij - (10) 

W 5 (W*. Wy ) * i 4 -ij 

and 


du 3 dx j . dv - il .v, , dw = dx 4 

Then transforming (6) from space variables to the lug variables, u. y. and w we have 


m ■ «sfffh , Ki(u+ X| Ki( u +v + x, ) (i(u + w + x, ) 
X 0(2cos0k)exp |-i s«i(i |k (u x + v x J - Akv x + kw x |J dudvdwdx, 


(II) 


Expanding the quadratic form K tJ k , k, in 0 and expressing the various lags in terms of 

u, v. and w, we have 

0 = exp^-4cos J 0 ^k 2 |2o J - K(u) - R(v) - R(u+v + w) + R(u + w) 

+ R(v-fvx) - R(w)| - kAk| 2o 2 - 2R(v) 

- R(u + v + w) + R(u + w) + R(v + w) - R(w)| 

(12) 

- kK | Jo- - R(u) - R(y) (■ R(u + w) + R(v + w) - JR(w)| 

+ k A k|o- - R(v) + R(y + w) - R(w)| 

+ (Ak) : |o- - R(y)| + fc~ |o- - R(wi| 
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On account of the moderate bandwidth assumption the leading term in the carrier wave- 
number k dominates the behavior ol the integral. Because of the deep phase modulation. ko»l, 
<t> is concentrated along two hyperlines in q. v, w hyperspace where the factor of k J is identically 
zero, f rom (12) it is seen that these lines are 


Cl. u = v * Q (w axis) 

C2. y. * - v, w = Q 


1 1 3 ) 


The thickness of the line masses, the e-tolding coherency distance, is (e.g., on Cl) 

lu c l 3 I v c I - (2kcos0o i )' 1 

where is the rms surface slope. Because <t> is concentrated about the two lines, we can de- 
compose M rito two separate mtegrals, M, and M 2 . about the small volumes surrounding Cl and 
C2. In each integral, the neighborhood of the origin, where the two line masses merge, will have 
to be ex luded. Consider th- integration about Cl, the w axis. Because of the dominance of the 
k- term, the cross-terms in k k, etc., and the (Ak) 2 term are negligible in the u, v integration: only 
the k 2 term that depends only on w will be important. Now. we can treat the problem in two 
ways at this point: we can allow for diffraction by the small-wave structure by assuming a com- 
posite surface consisting of two independent sets of smaller (diffracting) waves a.id larger (reflect- 
ing) waves. Analysis i then proceed by expanding <t> to first order in the small wave covariance 
function in the manner of Jackson (1974). The analysis, similar to that employed by Ruck, et aL, 
(1972) is complicated enough in two dimensions and it will be worse in three. To keep the three- 
dimensional analysis simple enough we will sacrifice what physical optics has to tell about the 
diffracted fields; accordingly we will stay strictly in the specular hackscatter regime where specular 
retro-reflection is the dominant mode of backscatter. (cf. Barrick. 1974; Brown, 1978) Thus, we 
will ignore the small structure and assume that the covariance function is well-behaved near the 
origin, i.e., that the surface is smooth on the order of a few wavelengths. Then the usual high 
frequency approximation can be made by representing the height covariance function R by its 
second-order Taylor series expansion about the origin (or about an^ lag w). 


Referring to (I’t, consider the VI integration over u ami v. On expanding the k’ term in u 
.iml \ about a fixed w t 0 to second order, there results the quadratic form, 

I » 

- — 2^ m ii , * >, i , j (141 

xv here the characteristic vector 

t (tj) 2kco*0(u x , Uy . v x . Vy ) (15) 

and the covariance matrix m l( ol the surface slopes. 



i^ 

1 

1 

1 

-J- 

1 

1 

-R xx (\y) - 

R x> (xv) 

m (J (w) 

1 

T| 

1 

1 

>. I 

r | 

i 

~R X y(w) - 

R nv (w) 


s\ mmetric 

°i 

°xy 



1 


°y 


where subscripts stand for partial differentiation, and 

-R xx (w) s <f x (0)f x (w)> (17) 


is the x-slope component covariance function. and similarly for \\ and yy. I lie slope variances 
(understood to be the radar-effective, diffraction limited slope variances, i.e.. having values less 
than the true high-frequency, optical values), 

o* -R XX (Q) h OS) 


With the s’ term' in xx that survives in the u, \ integration about VI xve then have for 0 

0 = c T . c -4cck 2 0k 2 |o 2 -K(w>| 


(l‘») 


Now, since* the coherency distances in u and \ are small compared to the extent of illumi- 
nation wc can set ja v Oin the gain pattern thus, (i, ... (I 4 - tc(\, X • ’ (x , 4 xx ) Also, m 
the pliasor in (I I) we can set Ak 0 since bx the moderate bandwidth assumption the phase dif- 
ference Aklv x l c «l. Iluis. if xve deline the joint specular slope vector 


N 


tanfld.O. 1.0) 


i '()) 


we can write the plnsor a\ c\p(-it • shrxp(-i2sin0xw x ). From the above, if we now integrate over 
the remaining space variable x, (which only appears in the gain function, aiul if we change from u. 
variables to t lie T variables we have for the integration over VI : 


M,(k) 


f.(w )e * co,2<,k2 I° J ' k ( 


. — — 1 

( 4 K 1 cos 4 0 ) J J. m 

{ ( ,r////e ^ .dt 4 |e*‘ J,mfl '‘ w xdw 


(21 > 


where 


■N 

i.l 


a ml A 


•/-uo 


I /*- 

w) - 1 (i*’l2i| K;- (x, ♦ 

\ J- uo 


( 22 ) 


U > d *l 


Note that the k-depcndence of M has disappeared in the high frequency k -*®» approximation. 

Tin* sign means that a region about the origin must be excluded in the integration ti e., where 

' I ^ 

VI ami V2 come together). We recognize that exp(- — } in^jtit, > is characteristic function ol the 

jointly, normally distributed slope vector £ (j' Xi , ( . f x , . fy, ) at two points I and 2 separated 
by the lag w I lie tenn in biaces, the Fourier transform of the characteristic function evaluated at 
the joint specular slope condition £ s' is the joint probability density function (pdf) ot ( evaluated 
al^| mz. (Fix/. I j». |frf)) 


Pls.yy 


\ = 


-Ini,,- Vi 


(2n)-’v/detm,j 


(2d) 


where m,, 1 is the inverse to the covariance matrix nij, given by (lb). 


Initiation (21) is written in the radar system coordinates with the x-axis lying in the plane 
of incidence. Define a new set of y coordinates such that the radar can assume any azimuth 


•1* relative to the new reference \-a\is (big. It. In (lux coordinate system, the surface modulation 
wavenumber 2ksiii0 becomes the vectoi modulation wavenumber, 

k B 2«t sui0( cos 4>. sin 4*); k s I k I (24) 

and the specular slope vector * becomes 

? * tan0(cos4*. sin4>, cos 4*. sin<t>) (25) 


Since the pdf (2.)) is invariant with respect to the orientation of the reference coordinate 
axes ( 21 ) becomes m terms of k 


M, (K) 


\. sec* 0 
IbR 4 



* cut * 0 1 u ^ - R (w»|. 


p(s.w V- ' -dw 


(2b) 


where now w = (w x , w\ ) is the lag vectoi in the fixed surface x. y coordinates. 


Unis far, vc have obtained, in the high frequency limit, the contribution M, to the moment 
M from VI. Since Mik. u. Ak) is symmetric with respect to a h and Ak interchange. and since M, 
is independent of Ak. and since M|(Aki 'M,(u) for arbitrary Ak * k it follows tit; ' 'he contribu- 
tion M, from V2. M. M,(k. Ak). rhus. the total moment M is given by (in the high frequency 
limit ): 

M(k. n. Ak) M | ( k ) + M, (Ak) (27) 


SI ( ONIM)RDI R I XPANSION Ol M, 

If x<rcot0 is not too large, the joint surface height characteristic function can be expanded in 
its argument, accurately, to second-ordei in the height covariance function R(w). The relevant 
parameter m this expansion is the significant wave slope, 

5 0 H s /l x , = (2/ir)K t> o (2S) 

where II,, = 4o is the significant wave height and l 0 is the wavelength ol the dominant wave cor- 
responding to the wavenumhei k„ of the peak ol the wave spectrum. In equilibrium sea condi- 
tions. " 1/20 - 5 , in strongly wind-driven developing seas the wave steepness may be as high 

N 

as 10' . (kinsman. l‘)bS). It we define the non-dimensional w avenumbei k k k 0 . then 
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( 2 «» 


ft = K <>cot0 - (v/2)A o Kcot0 

With values, /> 0 = 0.05. 0 - 15° and k -- 3. then 6 <0.88 and the error in the second-order ex- 
pansion will he less than 3 percent. 

It' the large-wave steepness 6 0 is small compared to the total radar-effective mis slops o,, 
then the off-diagonal Nocks of the slope covariance matrix m,j(w) for lags waway from the 
origin will he small compared to the diagonal blocks of m^. Then provided we are not too far 
into the wings ol the slope distribution, ie„ provided that tan 0/o, is not too large, we can ex- 
pand pts.wi about its value at infinite lag. where the slope correlations vanish, i.e., we can ex- 
pand about the value 

p(s;«>) s pf(s, ) (30) 

where the specular slope "sj --- tan0(eos<l», sin'M and p, is the (Gaussian) slope pdf (Cox and Munk, 
I ‘>54). The relevant parameter in the expansion of p is 

land 

e = 6 t , (31) 

A crude fitting of the geometrical optics backscatter cross-section. 

o° = rrsec*0P| ("sj ) (32) 

to the 2cm data ot Jones, et al.. (I ‘>77) gives the following estimate of the radar-effective rms 
slope at 2cm as a function of wind speed, W 

o j - 0.0025Wlms-‘| + 0.01 (33) 

This 2cm radar-etfective slope variance is approximately 60% of the total, or optical, slope var- 
iance ot Cox and Munk (l‘>54) over the wind speed range 10-25 ms' 1 . This is in perfect agree- 
ment with slope variances inferred from passive radiometric data ( Wilheit. 1979). With values, 
b 0 0.05,0= l5°,andW> 5ms' 1 , then i <0. lb 

l-oi typical nadir angles and sea conditions the parameters 6 and e are small enough that a 
linear approximation to M, should be fairly accurate, and that a second-order approximation to 


M| should bo very accurate. To lacilitatc the expansion of Mj we will assume that the x .uul y 
reference axes coincide with the symmetry axes of p, (s’, I. Then, if 

tantf 


cos>l' = a 0 cos4> 


land 

0 = sin 4* a 0 o sin4» 

On 


(34) 


p, ( s) becomes simply 


Pi * P|(«.0) 


1 t --jt° 2+ 0 2 ) 


;no x o y 


(35) 


Denote the waveheight directional spectrum F(K) l(K.*l>i. 

F(K) = ( 2ar)" 1 f K(w)c“‘— ‘-dw 
The surface height variance o' is the area under F(K), 

• lit 


(36) 


f Z n /* o© 

u- = / / F(k.<l')kJkd *1* 

A •'0 


(37) 


The directional slope spectrum k‘ F(K) is the contribution to the total slope variance o* front 
component waves travelling in various directions. 

•1 It />oo 

k 2 F(k.‘lMkdkd4> (3K) 


,• n 

•'ll 


k*F(K) is simply the Fourier Transform of the slope covariance function in the direction of 
analysis, i.e., in the direction of the component wave. Thus if the x-axis is again the tadar x- 
axis lying at an a/imuth relative to the fixed \-axis. we have 


K 2 F(K,4>) 


= (2tr)*‘ /- R 


xx (w x )c* ,kw *dw 


( 39 ) 


In the expansion of M ( to second order in the covariance functions and the consequent 
Fourier transformation in the surface reference coordinates we encounter such as 


/■ 


K x x e’ 1 * ' * d w - ( 2 * ) 2 cos 2 4>k 2 1 ( K ) 
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which follow# from the definition (36); the second order terms will result in various convolutions 
of (lie directional spectrum such as 



e-'K *Jw - C») J F(K)o K*F(K) 


(3*») 


where * denotes the two-dimensional convolution and the overbar notation will be used to denote 
groups of variables involved in the convolution operation. 


One last simplification will he made in the expansion of M, : we let the antenna beam in 
both range and azimuth dimensions be broad compared to the decorrelation distances of the sea 

A A 

;n either direction. I his allows us to set G(w) ~ 0(0) in the Fourier transformation over the 
modulation wave number range, and to set the rest of the integrand at its IK' value, p 2 ?s) in the 
Fourier transformation near k = 0. Further we will assume a separable pattern function G * 
Ci x tx)liy(y); then we have the definitions 


Lx 




‘txidx 


and 


Lx 


/ 


CiJ(x)dx 


(x - y) 


(40) 


A A 

From the definition (’2) u. the effective area A* we have; A e * LxLy. For Gaussian antenna 

Ac f 

beams, we have Lx = l.x/yj(x -» y) For simplicity we represent the Foruier transform of the an- 
tenna pattern simply as a DC' spike, i.e.. 



,.,-ikx 


dxdy -* 6(k) 


(41) 


From the above definitions, with the o 2 ,, elements set to zero in (Fq. (16)) on expanding the 
integrand ol 1 2M to secoiul order and Fourier transforming we have the result 


(\fl* LxLy \ , f 


v/7jt 

cot 2 fl( I - fi 2 ) 

Ly 


1 1 + (I + A) (a 2 + 0 2 ) 2 )K 2 F(K) 


cot 2 0 | K ’ F + + 2 a 0 K x Ky F + 0 2 K 2 fJ * FK 2 

^ ^ ^or 2 ( I - a 2 )k 2 1 * k- F + |4o 2 ar 2 0 2 k 2 F + a 2 ^ 2 ( I - (a - (1 ) : ) 

K^F + 4d‘ a 2 li 2 kfT | * kTk^F + 0 2 < I - tf 2 >k*F * k^F J } 


(42) 
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In (Ins we have let A - 6 4 /2(l - A 2 ). As mentioned earlier, the convolution (*) only involves 
those groups of variables carrying the overbar. 


Some features ol the solution should be noted. The most important is the reciprocal Ly 
dependence of the signal part (non-l)C’ part) of M,. The decrease in modulation signal power 
with increasing Ly is the price paid tor high directional resolution: of all the Fourier contrast 
waves contributing to the total modulation power the spectrometer is isolating a small subset that 
are travelling in the direction of radar look; consequently only a small fraction of the total spec- 
tral modulation power is being detected. The Lv'* dependence of tl.e solution can also be inter- 
preted as follows: In a real, short-crested sea. waves are running out of step with each other 
across the extent Ly so that modulation power will add incoherently as\/Ty. Thus the modula- 
tion power relative to the total average power decreases as L/ 1 *. and the relative spectrum of 
modulation power as Ly -1 . 


The directional resolution of the radar spectrometer is determined by ihe larger of the two 
effects (a) wave-front curvature on the surface and (b) the finite width of the data window 
Ly. The spectral window that results as the combined effect of (a) and (b) can be determined 

A 

by redoing the foregoing analysis in the Fresnel approximation. The lag vindow G(w) becomes 
in the Fresnel approximation (assuming a Gaussian gain pattern). 


G(w) 


G(wl 


exp 


/ l.x 2 w* + Ly 2 w 2 \ 
V 2irR 2 ) 


The corresponding spectral window ha>- a half power width. 


AK X = 


V3Rn: 


Lx 


I + 


k 2 Lx 4 esc 2 0~j 

CffR) 2 " 


1/2 


(x - y) 


The directional resolution is given by 


A‘h = 


Ak. 


(AKy « k ) 


(43) 


(44) 


I'or example, if R = 700km. k = 2*/ 200m. 0=15° and Ly - 10km. then 

A‘l* = 0.76(1 + 7.6) * 6.5 degrees. 


18 


I 


We note that the second order pjrt consists essentially of four convolutions involving the 
two-dimensional, directional wave spectrum F(K). In principle, since the second-order terms are 
small, we ought to he able to retrieve I from M, by a process of iterative deconvolution, start- 
ing with a first guess F(K) <* M, (K). 

To indicate, in rough terms, the amount of harmonic distortion resulting from the second- 
order terms in (42). vsc carry out a sample calculation in two-dimensions , using as a model spec* 
trum a one-dimensional Phillips* equilibrium spectrum, 

F(K) 3 BH(K - I )£* } (45) 

A 

where K = k/K 0 is the scaled wavenumber, and H is the lleavyside step function. We use a value 
of the spectral constant B = 0.01 as this gives a realistic value of the significant wave slope 6 0 
and hence, of the expansion parameters 6 and e. From the definition (28) of 6 0 and from (45) 
we have 6 0 = rr'v^B - 0.045. 

Figure 1 shows the results of the calculations for Q - 10° and 0=15° for a range of rough- 
ness conditions, i.e., for a range of the parameter a 0 = tanC/u, where o s is given by Equation (.13). 
For U - 10-15°, the harmonic distortion is quite tolerable, less than 30%, over the range of wind- 
speeds of practical importance. 


MEASUREMENT SIGNAL-TO-NOISE RATIO 

The result for M,. Equation (42), may be expressed as 

M,(K) = M,(0)[6(K) + P mod (K)| 

where l* mod is the modulation spectrum. 


(46) 


M, (0) is proportional to the square of the average backscattered power, we can express it 


as 


>/2* 

M. (0) -- | <IS(k H 2 >| 2 

Lx 


(47) 


The modulation index, defined for any a/imuth ‘l» is 

P(«l«) = f P mod (K. <*>)dK 


(48) 


l‘) 


NORMALIZED SPECTRUM NORMALIZED SPECTRUM 




0 12 3 4 


NORMALIZED SURFACE WAVENUMBER K 

Figure 3b. Kailar Spectrum vs. Surface Slope Spectrum 
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By (27). the total moment in 


M(k.c.Ak) - M,(0)|6(k> ♦ 5(AK) ♦ P mod (K) ♦ P mod (AK)) (49) 

\pm, we have k - 2xsmd (cos«I\ sin«l*». x oj/c, and similarly for AK. 


Analysis of Short l\ilso System 

The pulse bandwidth is assumed to be large compared to the dominant-wave modulation fre- 
quencies. Applying (4‘M to hquation (?). we find for the "per-pulse” (individual pulse) ensemble 
average output ol the spectrum analyzer. 


fi 2 

<Q(cj,)> = U ai — 
T 


I + p(4») 

r~ir 


+ P tnod(w 


. 4 ») 


(50) 


where the filter constant is set equal to I. 0 ai is the noise bandwidth of the H, filter; £ is the 
average energy of the backseat ter pulse. 


& 


2n <IS(k)l 2 > 



fc 0 (»»)l 2 dv; 


f is the effective backseat tered pulse duration. 


(51) 


T = L x/c 


(52) 


P c is the effective puls*, bandwidth, 


ft* - £>(w) - 


^F t ,(r)l 2 lK 0 (n - cj)l*dn 

J\ lE Q (i')l 2 dt' 2 


(53) 


Let r c = 2n, e/J k .(0) be the effective spatial pulse width. Then we can express the measurement 
Mgnal-to-noisc ratio, the ratio of the modulation spectrum to the broadband “Rayleigh clutter” 
spectrum, as 


4rrsin0 

SNR = P mod (K) 

T,. 


(54) 


where we have assumed u« I. Let the directional wave height spectrum by given by the 
Phillips* form 


I <k» = 5 x 10* 3 (8/3*) cos 4 4> k* 4 


(55) 


Let 0 • 1 5°; assume a nominal wind speed of 15ms* 1 and a corresponding dominant wavelength 
of 200 meters. With given h> (33), and assuming an upwind look, and using the linear part of 
(42) toi l’ lllo j we have P ni0( j * ,78m/Ly(km| and 

2.6 

SNR (56) 

r c |ni| I y|km| 

Reasonable SNR’s thus are possible on a per-pulsc spectrum analysis basis (i.e., without Doppler 
filtering) with r c on (lie order of a few meters and Ly on the order of 5-1 0km. SNR’s as low 
as -lOdB may be tolerable, if the measured spectrum is statistically stable, that is. if enough in- 
dependent pulses are averaged so as to produce a stable clutter noise spectrum 0J 1 relative to the 
modulation signal spectrum. The time-bandwidth product ( IBP) of the measurement is given by 
the analysis bandwidth )J a times the total integration time NpT where Np is a number of independ- 
ent pulses. At the c >0 f V confidence level, assuming large IBP. the length of the confidence inter- 
val is, relative to the signal level. 

l.64(SNRf‘ 

± — = £ KK’< confidence interval 

vfe l'Np/2* 

If we set the analysis bandwidth at 25 '. then, for our 200m wave example we have /J a T/2ir = 

(0. 25 K Lx/200 m). If SNR *0.1, I \ = 25km. and Np - 1000, then P mo j is measurable to ±9V? 
at the 90'. ' confidence level. 

As Alpers and llasseltnan (l‘>78) have shown for the narrow-band two-frequency measure- 
ment. Doppler filtering can be employed to increase the measurement SNR. One way to go about 
Doppler tittering in the short pulse spectrometer implementation is, in the time-domain, simply 
to keep track of or stare at the advected. frozen modulation pattern (e.g.. as in image motion 
compensation). The maximum integration time is determined by tr.e lesser of the wave period or 
the time it takes to move an appreciable fraction of the beam-spot in the along-satellite-track 
dimension. Both of these tailors will restrict stare times to something on the order of 1 second. 

If we allow a 50 < movement of the beam spot on the surface (in the along-track direction), then 
the maximum integration time available will be approximately, for looks to the side (4»> 30° ): 


f 


K A$cox<t> 



( 57 ) 


when* V' is tlu* spacecraft velocity. The Doppler spectrum width of backscatter Iron any fixed 
( stared- at i resolution cell on the surface is determined by the fastest interfering scattered waves 
Irom the extremities of the resolutk n cell ti.e.. at the 3dB azimuth gain points): 

V5v 

Alj„p sin0lsin4»l A4> (58) 


I lie last-time clutter spectrum dj 1 " ill be reduced by slow- tune Doppler filtering by the factor 
v/A i'j 0 pT, nl . Putting Ly - RA4\ we get for the measurement SNR after Doppler filtering, (again 
with 0 « 15°. •!»> 30° ) 


SNR £ 


r c v^R|km| X|m| 


a lesult which is independent ol the a/imuth beamwidth. With R = 700km. K = 0.02m 

7.4 

SNR * (60) 

r c | m | 


l ; or a 5 km l y. this is a !5dB improvement over the per-pulse SNR. 


Analysis ol I'w^v-i 'r equeue y System 

The transmitted waveform may be a sho't. wideband pitise or some tailored two-frequency 
wave! or m. Assume that the transmitted waveform Fourier transform t„(»') is constant over the 
bandpass d| and 0; ol the two filters ll ( and II, , then on applying (44) to Equation (3) we get 
for ensemble average output of the correlator, 


<Q i2 > = 


l\iiod(A , ’l 2 .4>) 


Die analysis bandwidth d a 1S the combined noise bandwidth of the two channels. If the H, and 
II, liber b. ndwidths p, ds then d a = • & •** the received energy in each channel: and T. 

again, is the effective back scattered pulse duration. 
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In (6 1 1, we note hut there in no background. Rayleigh “clutter” noise spec trum as there is 
in the short pulse system. This is because the expected value (or infinite ensemble average) of the 
cross-correlation of the basic Rayleigh clutter noise at two frequencies *»| and v } is zero (provided 
A»', J T»I). This follows front the definition of the basic unmodulated “clutter” process as a 
sample (of length T) of a stationary. Gaussian random process. 

The measurement noise in the two-frequency system is the residual clutter noise that sur- 
vives the cross-correlation over a l imte integr ation time T. If w c (t) is the (unmodulated) noise 
process going into the integrator (Fig. 2). then we have 

<0.) ■ / <w, t.t)>dt * 0 

•Ml 

Let a 2 be the variance of w c (t) and 0 C the bandwidth of w c (t). Then for large time-bandwidth 
products, 0 C T» 2rr. the standard deviation of the output ol the integrator will be Oq ^ 

y/2il0^f • o c T. 

A rigorous calculation of <Iq is beyond the scope of the present work; unfortunately. 

Patzen and Shiren (l l >5(>). in their analysis, stop short of calculating Oq for the two-frequency 
system. To get on with the analysis we will simply assume that o c = ti> 2 IT 2 and that 0 C * d a . 
Then, the per- pulse' SNR for the two-lrequency system will be on the order of 

SNR - (2»r ,/J i 212 T 1 ' 2 P mod (A»> 12 .4>> (62) 

The analysis bandwidth 0 a is set at some fraction of the modulation frequency in accord 
with desired spectral resolution. Let r a = (3 a /An, 2 . Then, using the same sample condition 
used in the short-pulse analysis Le.. wherein P mo( j(K.O) = 0.78m/Ly|km|, we have 

SNR ~ 0.055 r* /2 Lx ,/2 l km | Ly*‘ (km* 1 ! (63) 

Lett ini; Lx = 25km. ly = 5km, r a = 0.25. then SNR ~ 0.007. This figure is not as good as the 
per pulse SNR for the short pulse system. kg.. lor r c 2ni and I y >km, the short-pulse 
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SNK » 0.26, a factor of 32 greater. If the analysis bandwidth r is opened wide, r a * I, the two- 
frequency SNK becomes 0.06, and is comparable to the short-pulse SNR. The two-frequency 
SNR can also be improved by coheient Doppler integration. This would necessarily involve some 
kind of phase detection. If we assume an improvement factor ofVMdopTint (again with R ■ 
700kn», X = 2cm), then the measurement SNR will be 0.1 (r, a 0.25) and 0.85 for r, * I. (The 
reader is cautioned that we have been loose with certain numbers, i.e., with the effective integra- 
tion time, Tj n( and Af<j 0 p, and o c .) 

( ONCLUSION 

We have presented a three-dimensional physical optics solutior for the generalized fourth- 
order moment of the surface transfer function in the case of near-vertical, specular backscatter 
from a Gaussian sea. and we have applied this solution to the analysis of the short-pulse and 
two-frequency techniques for the remote measurement of ocean wave directional spectra. The 
predicted modulation spectrrm is found to bear a fairly good fidelity to the directional slope 
spectrum; harmonic distc ■' n, or spectral smearing, is on the order of 20-30% or less. The 
solution’s only serious defec , it is felt, is the Gaussian assumption for the ocean wave statistics. 
The measurement signal-to-noise problem for both short-pulse and two-frequency spectrometer 
implementations was examined briefly, and it was found that satellite measurement of the dr- 
.-ctional wave spectrum with either of these spectrometer systems appears to be feasible; it remains 
for future work to model diffraction and non-Gaussian wave statistics, and to perform a more 
thorough measurement signal-to-noise analysis; for example, trade-offs involving transmitter power, 
antenna gain and thermal noise must be considered. 
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